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ABSTRACT 

i'lxed configuration feedback control often forms the 
basis for the design of a linear control system. The 
design freedom is limited to the adjustment of the free 
parameters of the system. Analytical methods prove to be 
valuable in these circumstances and the aim is the optimi- 
zation of a selected index of performance. A suggested 
index is developed for a regulator system from a comparison 
of the dynamics of the subject sy4ftem with those of a 
desired model along a transient state space trajectory. 

The result is an optima,! system dependent on the initial 
state chosen for the transient response. The effects of 
the initial state on the optimal system are Investigated 
and the final selection is based on statistical consider- 



ations 
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CHAPTER I 



INTRODUCTION 

The design of feedback control systems sometimes in- 
volves the selection of free parameters in the plant and/or 
compensator in order to satisfy the requirements of the per- 
formance specifications. The configuration of the system is 
assumed to be fixed, i.e., a compensator has already been 
chosen and the design freedom is limited to the adjustment 
of one or more parameters. The trial and error process for 
the selection of several parameters represents a very labor- 
ious and sometimes impossible task for the classical design 
by root locus methods. Analytical methods are appropriate 
in this case. A measure of quality is established for the 
system in question by the construction of an index of per- 
formance, xvhich upon minimization (or maximization) will 
yield the optimum values for the free parameters of the 
system. 

Different criteria could be adopted for the construc- 
tion of tae index of performance, all normally based on 
physical reasons and each satisfying particular requirements. 
Many performance' indices originating from good engineering 
reasoning lead to analytically insolvable problems. A com- 
promise in the selection of the index is then necessary to 
obtain a solution. The classical forms of the performance 
specifications, such as rise time, maximum overshoot, 
settling time, etc., in the time domain or bandwidth, phase 
and gain margins, etc., in the frequency domain, are trans- 
formed into a single measure by conceiving an ideal model 
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whose behavior meets all the performance specifications, A 
mathematical relationship established between some variables 
of this model and the corresponding variables in the actual 
system serves as a measure of quality as the free parameters 
are changed in a search for the optimum adjustment. The 
time integral of the difference-squared between the output 
of the model and the output of an actual system is a practi- 
cal index of performance when the system is subjected to 
transient input signals. Here the model has a transfer 
function equal to unity. This is usually solved with the 
aid of Parseval*s Theorem [23. 

An index of performance has been derived for the regu- 
lator problem[^3. This index of performance, or cost func- 
tion, is a quadratic function of the state vector. 




( 1 ) 



The matrix of the quadratic form, Q, is the weighting 
matrix and contains the free and fixed parameters of the 
actual system and some weighting factors determined by the 
model or ideal system. The model obtained from the per- 
formance specifications is of the same order as the actual 
system. The weighting factors are adjusted in such a way 
that if all the parameters of the actual system were free, 
the cost faction would exhibit an absolute minimum at the 
model values. If the actual system is constrained to have 
less than n independent free parameters, a minimization of 
the cost function will yield the optimum values for these 
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parameters. 

The measure of quality is obtained by disturbing, the 
system with some initial conditions. For a stable system, 
the state space trajectory ends up at the origin for large 
values of time. This type of excitation is characteristic 
of the regulator problem. By minimizing this cost function 
the trajectory of the subject system is forced to approximate 
as closely as possible the trajectory of the model. If all 
n parameters of the system were free, it would Indeed simu- 
late Identically the dynamics of the model and hence satisfy 
in full agreement the performance specifications. In any 
other case, the model behavior will not be reproduced, but 
the actual system will attain its best performance as 
measured by this optimization criterion. 

The importance of this method over the common types 
of indices is the consideration given to every part of the 
dynamics of the system and not to a single manifestation as, 
for example, the error in position. The contributions from 
all the state variables are properly weighted according to 
the characteristics of the model, which carries the infor- 
mation of the performance specifications. 

The evaluation of the Integral in equation (1) is done 
by Liapunov functions with vector-matrix methods, and the 
final form is shown to be a quadratic function of the 
initial state, 

J = x^(0)Px(0) (2) 

The matrix of this quadratic form, P, is related to the 
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weighting matrix and the actual system’s matrix. 

The nature of this solution implies that the cost 
function will be different for different initial states and 
hence the optimum system will depend on the initial condi- 
tions chosen for the calculation. 

In this paper, a similar index of performance is de- 
veloped, based upon a different approach than that described 
above. To obtain a measure of performance, at every point 
on the state space trajectory of the actual system, a com- 
parison is made between the dynamics of this system and the 
dynamics the model would have if it were at that particular 
point. 

The index of performance can be expressed in the same 
form as equation (1) and hence is also a function of the 
initial state. However, the weighting matrix is now direct- 
ly obtained from the coefficients of the model and actual 
system characteristic equations. This is a much more direct 
method than that for obtaining the Q matrix in equation (1), 
which relies on complicated algorithms for the weighting 
factors. In the analysis of this index it is found that 
the results are exactly the same as for the index of equation 
(1), Further Investigation shows that both cost functions 
are intimately related by a constant factor and an additive 
constant. The effects on the optimum system of the initial 
state chosen for the calculation are investigated in this 
paper and the final decision on the design is based on the 
statistical characteristics of the disturbing initial con- 
ditions. 
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The Hex'; concept developed also suggests a time-varying 
optimization scheme by which the actual system would be made 
to follow exactly the state space trajectory of the model by 
continuously varying its free parameters so as to match at 
every point the dynamics of both systems. Examples are pre- 
sented to complement the theory. 
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CHAPTER II 



THE INDEX OF PERFORMANCE 

The first step in the process of analytical design is 
to establish the concept on which the index of performance 
will be based. The objective of this performance index is 
to include in a single mmber a measure of quality for the 
performance of the system. Factors such as mathematical 
difficulties will perhaps force the designer to modify his 
performance index in favor of mathematical simplicity, hope- 
fully keeping the essence of the idea that originally moti- 
vated it. Once the configuration of the system iq fixed and 
the index of performance has been defined, the design 
proceeds in a straightforward manner. For the type of index 
to be considered here, the performance specifications are 
transformed into an ideal or model system thus fixing the 
location on the s-plane of the required zeros and poles of 
the transfer function. This model is of the same order as 
the actual system and fully satisfies the specifications. 

The model and actual systems can be represented by the signal 
flow graphs shown in Figure 1. 

The feedback coefficients of the model are defined by 
the performance specifications. Some of the coefficients 
of the actual system are constrained by the physical nature 
of the plant, and the others are functions of the free para- 
meters. If all coefficients in the actual system were free 
it would be trivial to reproduce the model. Generally this 
is not the case and the free parameters must be adjusted 
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according to some criterion, 

A measure of quality vrhich originates from a good 
physical concept is given by 



oo 




xirhere the integrand is the square of the difference between 
the n derivatives of the model and actual systems differ- 
ential equations respectively, when the system is excited by 
initial conditions only. This index gives a comparison of 
the complete dynamics of both systems in terms of the high- 
est derivatives, and does not require any weighting factors 
since it implies a direct relationship between the two 
systems. The minimization of equation ( 3 ) as function of 
the free parameters results in an actual system verbose dynamic 
behavior is a close approximation to that of the model. 
Hovrever, the development of this performance index into a 
practical form presents the impediment of requiring the 
solution of both system differential equations, which shadows 
the merit of its physical meaning and renders the integra- 
tion of equation ( 3 ) very difficult. It is then necessary 
to conceive a similar structure for the index of performance 
but one which permits convenient mathematical manipulations. 

A similar way of evaluating the performance of the 
system can be obtained from a comparison along a trajectory, 
of the n^'^ derivative of the actual system differential 
equation with that of the model when both systems are in the 
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same state. It Is again assumed that the system is driven 
by initial conditions only. The motivation for this index 
arises from a consideration of a second order process. Con- 
sider the system shown in Figure 2. The system differential 
equation is 

X + a2X + aj^x =0 (Zf) 

where a^^ is the free parameter 

The performance specifications require a model system char- 
acterized by the differential equation 

y -I- + Siy =0 (5) 

The actual system trajectory in the phase plane is shox'm in 
Figure 3* At any point on this path, x is determined by the 
coordinates of the phase plane and the coefficients of the 
differential equation 

x^ = (^>) 

Also, equation (5) determines the value of y at any point on 
the phase plane 

^ "^ 1^1 

vihere the subscript 1 Indicates the particular value of the 
variables at any given point. 

Along any trajectory of the actual system there will be 
a continuous function of the state giving the value of y. 
This function is 
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model 



actual system 



Fig. 1 • Signal flow graphs for model and actual system 



i 




Pig. 2. A typical second order system 



1-5 



y = 



( 8 ) 



- SjX 



The optimization criterion is the minimization of the differ- 



xmen actual system and model are in the same state. Then 



Figure k shows a signal flow graph for the calculation of 
the index. 



trajectory along x^rhich this comparison is made will change 

until the value of the index is a minimum. This will be the 

closest trajectory to that of the model and the performance 

of the actual system will have been optimized according to 

this criterion. In terms of the dynamics of motion, we 

vrish to make the acceleration of the actual system at any 

point on the phase plane a best approximation to the accel-^ 

eration the ideal system would have at that point. 

t Vi 

This idea is extended to an n order system 



ence squared between the variables x and y along the actual 
system trajectory. The index of perforaance becomes 




( 9 ) 




( 10 ) 



& 



In searching for the optimum adjustment of the 




+ • • • + = 0 
at 



( 11 ) 
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Pig. 3. Phase plane trajectory for system of Pig. 2 




Pig. 4, Signal flow graph showing the structure of 
the index 
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where the index of performance becomes 




when actual system and model are in the same state. 







By defining the state vector 
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this index can be written in quadratic form 

J- =y [iya„fxl + + • • • 

“ [Xi - - - Xy,] 
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v/here bj^j = (Sj^ _ ) (Sj - ) , i = 1,2, 

n. 



, n and 



0 . 2 ) 

Then 

(13) 

(14) 



(15) 

(16) 

j = 1,2, 
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Then 




(17) 



vjhere the matrix of the quadratic form, B, is a positive 
semi-definite, symmetric matrix. 

This is the general form of the index of performance 
as given in Equation (1), It is solved hy vector-matrix 
methods using Liapunov's Theoremt^J. The solution is 



J = x^(0) P x(0) 



( 18 ) 



where p is the matrix of the Liapunov function and is re- 
lated to the weighting matrix B and the actual system matrix, 
Aa, by 



-B = + PA 

a a 



( 19 ) 



The B matrix is directly obtained from the model and 
actual system coefficients as shovm in equation (16), The 
actual system matrix is defined by the matrix differential 
equation in canonical phase variable form 



vmere 




. 0 ± 0 - -- 0 
0 0 1 

I 

I - ^ 

N. 

' V 1 

--- _aa 



( 20 ) 



(21) 
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Signal flow graphs are helpful in solving equation 
(19) for the elements of the p matrixL^J. The results for a 
second order system are stated here for illustration. 



R. 



- - ^1?, + 






Zdi, 



bii + b 



2da 






P-2 " ~ 



( 22 ) 



P = ~~ ^11 -h — ^11 

rzz zaz ^ Zdxcxi " 

The cost function, equation (18) becomes 



J — X ]_(0)Pj^] + 2Xj^ ( 0 )x^ ( 0 )pj^2 *t ^2 ^®)P22 (^3) 

The minimization of the cost function vrill yield the 
optimum adjustment for the free parameters. Because of the 
nature of this artifice, the optimum system will depend on 
the initial state chosen for the calculation. This problem 
is discussed in Chapter III. 



A TIME-VARYING OPTIMIZATION METHOD 

The minimization of the index of performance developed 
in the last section produces a set of parameters for which 
the dynamics of the system is a best approximation to that 
of the model. In terms of the state space, this means that 
the actual system trajectory is now the closest trajectpry 
to that of the model which can be achieved with this con- 
strained system under the effect of a given set of Initial 
conditions. The index is nov; a measure of the total error 
incurred by the system in going from this particular initial 
state to the origin of the state space. Becaiise it is a 
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quadratic function, it gives more Importance to large errors 
than to small errors but is not capable of reducing the mag- 
nitude of any particular error in performance. 

The criterion on which this index is based suggests a 
possible way to obtain with the actual system an exact re- 
production of the model dynamics if the free parameters are 
allowed to vary continuously according to some function of 
the present state. Actually, it would be enough to have 
only one free parameter in an n order system to implement 
this scheme. 

Consider again the second order system of Figure 2. 

The actual system and model differential equations are given 
in equations (4) and (5) respectively, vifhen both systems are 
in the same state. 

y = X 

• » 

y = X 

and the condition for the tvjo state space trajectories to 
be the same (and so the dynamics of both systems) is 

•« «• 

y = X 



or 



g^x + ggX = a^x + a^x 

then 

"I" ^ ^ 



(24) 
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Equation (24) is the required functional relation be- 
tween the free parameter and the state of the system. This 
is plotted in Figure 5 for a typical case. 

The curve shows a discontinuity whenever the trajec- 
tory crosses the x-axis in the phase plane. Since infinite 
gain is not practically feasable, the function can be approx 
imated as shown in Figure 5. A change in the sign of the 
gain factor from time to time is also required. The reason 
for the discontinuities is that when x = 0, the free param- 
eter, a^ , which modifies the variable x in the feedback path 
has no effect on the behavior of the system, now solely con- 
trolled by the variable i and its fixed feedback parameter 
a^. Thus, an infinite gain is the result of the impossi- 
bility of altering the system's trajectory. 

th 

For the general case of an n order system, the con- 
dition under vxhlch the actual system will duplicate the 
model dynamics is 






X + 



h 



c/x 

c^t 



n, 









(25) 



If we assume that a^ is the only one free parameter, then 



Equation (26) shows that only one free parameter is 
needed to reproduce the model dynamics if the state of the 
system is known. Normally the fact that the state vector 
is completely measurable implies that the feedback paths 
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Pig. 5. A plot of k(x) versus t; model: y + 2.8y + 4y = 0 
actual system: x + x + kx = 0 and x^(0) = (1 0) 
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can be modified bo copy the model characteristics exactly. 



This discussion would apply in those special 
the state vector is measurable but cannot be 
in the feedback channel. 



cases in v:hlch 
used directly 
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-CHAPTEH III 



THE OPTIMUM SYSTEM AS A FUNCTION OF THE INITIAL STATE 

The index of performance that will yield the optimum 
system is a function of the initial state, Ne wish to in- 
vestigate the effect of the initial conditions on the 
optimum pole locations. These optimum poles will occupy 

certain regions on the s-plane according to the constraints 

; ) 

of the system and the requirements for stability. Also of 

i 

interest is the knowledge of hov/ the extension and shape of 
these regions vary as a function of the relative position 
of the model poles vilth respect to the actual system’s 
attainable root regions. 

Only the normalized initial state vector will be con- 
sidered in this discussion since all other points in the 
state space are proportional to this unit vector and differ 

in cost function only by a constant factor. For each 

^ > 

orientation of the unit Initial state vector there is a 
corresponding optimum set of poles for the actual system 
and a value of the cost function which is a minimum for 
this particular choice of initial conditions. Any other 
orientation will yield a different set of optimum poles. 
However, it will be seen that, because of the symmetry 
exhibited by the index of performance, only a half of the 
state space need to be searched for optimization. The 
design must provide a system whose performance is optimum 
in some modified sense for all possible initial condiijions. 
The selection of such a set of optimxim parameters must 
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be done on the basis of the statistical properties of the 
disturbance. If the probability density function can be 
established, then the average value of the cost incurred by 
the system can be found for all sets of optimum parameters. 
The set that provides the minimixm average cost is the 
design selection. This, of course, implies the consider- 
ation of the behavior of the control system as a regulator. 
It is not possible to obtain an analytical relationship for 
the optimum poles as a function of the initial state and 
the proximity of the desired model poles for a given con- 
strained actual system. However, the general trend is that 
the regions containing the optimum poles reduce, and finally 
collapse to a point, when the model poles get closer to and 
reach the permissible regions for the actual system. Simi- 
larly, the same effect is obtained by gradually releasing 
the constraints of the actual system but keeping the same 
characteristics of the desired model. These effects will 
be discussed for the second and third order systems by 
means of examples, 

SECOND ORDER SYSTEM 

The results of the optimization scheme developed in 
Chapter II will now be discussed for a typical second order 
system. Two parameters define the dynamics of this system. 
The method is applicable when one of this parameters, 
usually in the plant, is fixed. Nlth only one degree of 
freedom, the possible poles for this system must lie on a 
definite locus. A pair of required poles outside this locus 
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is simply unattainable and the optimization process must 
be employed. 

Consider the second order system shown in Figure 6. 
'The transfer function is 



■ C ^ ^ k (2?) 

^ + s + k 

The free parameter is the gain of the system TJhereas tne 
plant contains the fixed parameter. The performance 
specifications require a model defined by the transfer 
function 

ic , ;'v , 

C ^ 4 (28) 

' ^ 's^ + 2.8s + 4 

It is noted that the model poles are unattainable by 
the actual system. It is then necessary to find an adjust- 

t’ : • t 

ment for k that will give the best performance with respect 

to the model according to the criterion previously establish- 

! . . ■ • . , . .. , ■ 

ed. 

The cost function is given by equation (23) 

J = X£^(0)p^^ + 2 x^(0)x2(0)p^2 



The' elements of the p matrix are those given in equation 
(22). The weighting matrix from equation (l6) is 



(k-4) (k-4)(l-2.8) 



B = 



(l-2.8)(k-4) (1-2.8)' 



(29) 
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then, J becomes 



J = x^'^(0)(0.5k -0.08k-3.2-i-2) + 2Xj^(0)x2(0)(0.5k-4+|) + 



X2^(0)(0.5k-2.38+|) 



( 30 ) 



Only points on a unit circle in the initial condition 
plane need be Investigated. Equation (23) is symmetric with 
respect to a diameter on the initial condition unit circle. 

As a further simplification, only the upper semicircle will 
be considered. 

The root locus for this system is shown in Figure 7. 

For each value of 0, where 0 is the orientation angle of the 
state vector around the upper unit semicircle of the initial 
conditions plane. Figure 7(a), a value of is found by 

minimizing the cost function with respect to this parameter. 
The plot of optimum regions for the corresponding poles of 
the actual system is shown in Figure 7(t>) together with the 
desired model poles. The extremes are located at (-0.5+ 
j 1 . 936 ) and (-0.5+J 0.775) ^or the upper half of the s-plane. 
Each pair of conjugate poles on this region will give an 
optimum performance only when the system is excited by the 
initial conditions that yielded that pair of poles. For 
any other initial state there will be a different set of 
optimum poles and a correspondingly different setting of 
the free parameter for optimism response. It is not possible, 
of course, to modify the adjustment of the free parameter 
according to the type of disturbance or even worse, to pre- 
dict its values. It is then required to adopt a single 
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Fig* 6, A constrained second order system 




model 
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Fig. 7* (a) Initial condition plane, (b) Root locus 
for s2 + s + k = 0 
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value for the free parameter that will give the best average 
index for all possible situations. Figure 8 shows the value 
of the minimum cost as a function of initial conditions. 

The horizontal axis represents the orientation of the Initial 

state vector. This curve shows that a vector orientation 

i 

of 140*^ will give an™ absolute minimum for the cost function 
of this example. This corresponds to kQp^ = O .85 and the 
optimum pole lies right at the bottom of the allowable 
region. This would seem the best adjustment to use in the 
actual system. Figure 9 shows the cost incurred by the 
system as a function of the initial state for different 
optimum settings of the free parameter. It is seen that 
the optimum parameter obtained at the absolute minimum of 
the cost function of Figure 8 yields the highest index 
when used with all initial conditions. The area under 
each curve is a measure of the cost accumulated by the 
system when excited by all initial conditions on the upper 
semicircle of Figure 7(a) for a given optimum adjustment 
of the free parameter. Figure 10 shows a plot of the 
total cost as a function of the free parameter. 

The curve in Figure 8 also exhibits another local 
minimum, which for some models is the absolute minimum of 
the cost function. The corresponding optimum parameter 
also produces a high average index as shown in Figure 9* 

The minimum area is obtained for a = 2.39> which 

corresponds to an initial state vector orientation of 25 °. 
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Fig. 6. Minimvun cost versus initial conditions 




Pig. 9. Cost as a function of the initial state for diffe 
rent optimum settings of the free parameters 
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This particular value of k would be the design selection 
when there exists a uniform probability on the initial 
s ta t e 



p(0) = i de , 0 ^ 180° 

TT 

For any given probability distribution on the initial 
state, p(0), we define a modified index of performance, 

X 

which will be independent of the initial conditions, 

Jmod = / P(e) J(e»k) dO ( 31 ) 

o 

The minimization of this index with respect to the free 
parameter yields an adjustment for which the average cost 
of the system is a minimum. As a consequence, the maximum 
value of the cost for any possible initial state will also 
have been minimized. This is equivalent to a minimax 
solution in optimal control theory 

For the case of the uniform probability considered 
above, the modified index is 




The minimization of equation ( 32 ) with respect to k yields 
the result previously obtained. 

The size of the optimum region depends on the prox- 
imity of the model poles to the root locus of the actual 
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system. Figure 11 shows how these regions increase as the 
model poles get further apart from the root locus of the 
actual system. 

Only the upper half of the s-plane is shovm. For 
model poles on the root locus of the actual system the 
regions are simply reduced to a point. The growth of the 
regions as the model poles get further apart from the actual 
system root locus is more accentuated in the positive jco 
direction, an indication of the tendency of the index to 
make the natural frequency , caJ^ , of the actual system equal 
to that of the model. The effect of the value of the 
imaginary part of the model poles on the size of the 
region is not pronounced, except that for small values 
of cj , some optimum poles could lie on the real axis. Only 
near model poles are of interest here since remote poles 
would probably require a change in the configuration of 
the actual system, x^hlch is outside the scope of this 
discussion. 

The results stated so far are all based on the minimi- 
zation of a performance index which is a function of the 
free parameter for a given set of initial conditions and 
can be represented by a two-dimensional curve. This curve 
is actually the Intersection of a plane and a more general 
cost surface obtained when both parameters in the actual 
system are free If in equation (23) the actual system 

parameters are variable the cost function becomes 
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Pig. 10, Total cost as a function of tia© fre© paramator 




Pig. 11. Size of optimum regions for several models 
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(33) 



The shape of this cost surface will depend on the 
initial state used, as shown in equation (33) but the mlni- 
mum value of this cost will always be zero and lie on the 
point (4,2.8), which corresponds to the model of this 
example. The cost surface is shown in Figures 12 and 13 
for two different values of the initial state vector orien- 
tation. It can be seen that while the absolute minimum of 
the surface is fixed, the local minimiun of the inter- 
ception of the plane a^ = 1 with these surfaces changes 
with initial conditions. 

THIHD ORDER SYSTEM 

Consider a third order linear control process charac- 
terized by the differential equation 



eters to a value of 2. The specifications suggest a model 



X + 2x -»• a^x + a^x = 0 
2 1 



(3^) 
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12* Isometric projection of cost surface for a model 
characteristic equation s^ + 2.8s +4=0 and an 
initial state vector x^(0) = (1 0) 
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Fig. 13* Isometric projection, of cost surface for a model 
characteristic equation s^ + 2*8s +4=0 and an 
Initial state vector x^(0) = (0*866 0*5) 
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given by 



y 7.8y + l8y + 20y = 0 



(35) 



Because of the constraint in the actual system, these 
characteristics cannot be duplicated. The free parameters 
have to be adjusted to obtain the optimum performance. The 
cost function from equation (18) is 



J = 2x^(0)X2(0)p^ 2 2x^(0 )x^(0)p^^ + 

X2^(0)p 22 't- 2 x2(0)x^(0)p 2^ + x^^(0)p^^ (36) 



where p^j is a particular element of the matrix p obtained 
from equation (19)^^^. The weighting matrix is now 

2 



B = 



(a^-20)‘ 



(a2-l8)(a^-20) 

(2-7.8)(a^-20) 



(a^-20)(a2-l8) 

(a2~l8)^ 

(2-7.8)(a2~l8) 



(a^-20)(2-7.8) 

(a2-l8)(2-7.8) 

( 2 - 7 . 8 )^ 



(37) 



The cost function is then 



J = 3C| (oi 
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here D ; 


= ^2^3 " 


ai 












The 


effects 


of 


the 


init 


ial state 


on the optimum 


pole 



location snoild be investigated for all points on tne unit 
sphere in the state space. Because of the s/.nnetrj of 
equation (36) with respect to a diameter, only a hemisphere 
need be considered. For t'.fo degrees of freedom, the cost 
function corresponds to a three-dimensional cost surface. 

The shape of tnis subsurface nill change vith initial state 
and its corresponding minimum point v^ill move around on the 
plsne. The subsurface coi'respoading to a coefficient 
= 7. '3 .’/ill have .a minimum value of zero at the model 
values, independently of the initial state. A typical sub- 
surface is shov/n in Figure 14 for the case x^(0) = (-0.754 
0.133 0.643). By computer searching techniques the optimum 
coefficients for this system uere found for several initial 
conditions on tne unit he.iiisphere to the right of the x-z 
plane . 

The roots obtained from these optimum characteristic 
equa Cions conglomerate in limited regions on the s-plane. 
Tnese regions are sho-./n in Figure 15* Observe that the 
problem of stability has been automatically solved. If the 
Gonstrai.it on the actual system 7/ere released these opti- 
mura regions ./ould collapse into points './here the model poles 
are and the cost function would be identically zero for all 
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Fig* 14* Isometric projection of cost subsurface for a mo- 
del characteristic equation s5+7,8s^+l88+20=0 and 
an actual system with a constraint a^=2* The 1- 
nltlal state vector Is x^(0)=(-0.754 0.133 0,643) 
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0 Optimum poleB for 
x^(0) = (1 00) 

O Optimum poles for 
modified index 

X Poles by direct re 
production of model 
coefficients in sys 
tern free parameters 




Fig. 15. Optimum pole regions for third order system with 
constraint in = 2 
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initial conditions 



Each set of optimum poles is associated with a value 
of the cost function, which is a minimum for the initial 
state from which this set was obtained. For any other 
initial state there will be a different set of optimum 
parameters. The final selection of parameters is based on 
the statistical characteristics of the disturbance. The 
modified index is obtained from an extension of equation 
( 31 ) to a three-dimensional space. For a uniform proba- 
bility distribution over the hemisphere this index is 



This modified index is independent of the initial 
state. The minimum average cost function, and hence the 
minimax solution is obtained for this example with the 
optimum coefficients 



The corresponding poles are shown in Figure 15 . 

A reasonable choice is to assume an Initial state 
vector 



The optimum set of poles for this case is shown in Figure 





(39) 



= 6.5 

a£ = 8.7 



x'^(O) =(10 0) 
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15. The cost subsurface is shown in Figure l6. This set 
of initial conditions resembles the widely used step 
function in control systems design. 

A x^rorthwhile observation is that the available free- 
dom in the system parameters could be used to reproduce 
exactly some of the corresponding model coefficients. Con- 
sider the actual system and model differential equations 
as given in equations (3^) and (35) 

X + 2x + a^x a^x =0 (34) 

y + 7.8y + I8y + 20y =0 (35) 

The direct reproduction of the model coefficients in the 
actual system will yield 

X 2x + I8x + 20x = 0 (40) 

The corresponding poles are close to the optimum regions 
obtained by the previous method. They are shown in Figure 
15. This possibility in parameter selection is limited by 
the stability requirements a^>0, a2>0, a^a^-a^ >0. vi/hen- 
ever possible, this criterion gives a good, time-saving 
procedure if the third requirement for stability is well 
met. 

It is interesting to see how the optimum poles migrate 
as the Initial state vector sweeps the normalized state 
space. A mapping of a unit circle on the x^-x^ plane in 
the state space onto the s-plane is shown in Figure 1? for 
this example. The optimum poles return to the same point 
after the initial conditions vector turns 18O®. 
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Fig. 16. Isometric projection of cost subsurface. Model 

and actual system as in Figure 14. Initial state 
vector i^lO) = (1 0 0) 
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Pig* 17. Mapping of the unit circle on the 
to the 8>plane 



Consider nov; taat for tne same performance specifica- 
tions, the constraint on the actual system is somewhat less 
stringent. Assume that the system's cnaracteristic equation 
is 



The model poles will novj be closer to the attainable 
regions of the given system and consequently, the results of 
the optimization procedure will be Included in smaller re- 
gions. These regions are shown in Figure 18. There is also 
a net displacement of these regions to the left and now the 
complex conjugate model poles are included. The set of 
poles obtained by the mere duplication of model coefficients 
in the free terms is also shown. 

Using the modified index of equation (39) foi^ a uniform 
probability on the initial state, the final selection is 



s 




+ + a^ = 0 



(^ 1 ) 



a 



1 



16.3 



= 15.2 



The optimum poles are shown in Figure 18. 
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Pig. 18. Optimum pole regions for third order system 
with constraint in = 6 
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CHAPTER IV 



EVALUATION OF RESULTS 

Because of the constraints on the actual system and 
the requirements for stability, the poles of the system 
transfer function are restricted to certain regions on the 
s-plane. The optimum pole regions have a size, shape, and 
position according to the proximity of the model poles to 
the alloviable regions. The selection of system parameters 
for optimum design was based on the statistical properties 
of the disturbing initial conditions. That set of poles 
from the optimum regions which gives the minimum average 
index imder all possible initial conditions is the desired 
solution to the problem. 

The index of performance developed in Chapter II v/as 
given the form of equation (1) for mathematical convenience. 
The tv^o indices differ in the description of the weighting 
matrix. However, the optimum parameters obtained are 
identical for every initial state. The elements of the Q 
matrix in equation (1) were chosen to place the absolute 
minimum of the cost function at the coordinates of the 
coefficients of the model characteristic equationt^^"- while 
the elements of the B matrix in equation (1?) came directly 
from the general expression of the performance index of 
equation (12) when expressed in matrix form for mathematical 
manipulation. 

A comparison between the two indices is done for the 
second order system. The weighting matrices, and B, 
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The expanded form of the cost function given in equation 
(2) is, by use of equation (22) 
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Where = index of performance with matrix Q 

J = index of performance with matrix B 

B 






'*9 



Equations (44) and (45) are related by 

~ ^ ■*■ 2x^(0 )x^(0)g^ + x^^(0)g^ (46) 

For any given initial conditions, the right hand side of 
equation (46) is a constant related to the model coeffic- 
ients and the particular set of Initial conditions. Thus, 
in general, both cost functions differ by a constant factor, 
which is the square of the model coefficient g^, and an 
additive constant given by equation (46). As a conse- 
quence, the optimum coefficients of the actual system 
obtained by the minimization of these two Indices are 
Identical. This is true for systems of any order. How- 
ever, the algorlt'nm to determine the B matrix is simpler 
than that for the Q matrix. 
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